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Abstract 

Cosmic high energy neutrinos are inextricably linked to the origin of cosmic rays 
which is one of the major unresolved questions in astrophysics. In particular, the high- 
est energy cosmic rays observed possess macroscopic energies and their origin is likely 
to be associated with the most energetic processes in the Universe. Their existence 
triggered a flurry of theoretical explanations ranging from conventional shock accel- 
eration to particle physics beyond the Standard Model and processes taking place at 
the earliest moments of our Universe. Furthermore, many new experimental activities 
promise a considerable increase of statistics at the highest energies and a combination 
with 7— ray and neutrino astrophysics will put strong constraints on these theoretical 
models. The detection of ultra high energy neutrinos in particular is made likely by 
new experimental techniques and will open an important new channel. We give an 
overview over this quickly evolving field with special emphasize on new experimental 
ideas and possibilities for probing new physics beyond the electroweak scale. 
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1 Introduction 



The furtherst continuous source so far observed in neutrinos is the Sun ffl. These solar 
neutrinos together with their atmospheric cousins [[J, |3| which are produced by cosmic ray 
interactions in the atmosphere have in fact provided evidence for neutrino mass and oscilla- 
tions. They are thus playing an important role in probing new physics beyond the Standard 
Model. The only extragalactic source so far observed in neutrinos in the 10 MeV range was 
sporadic, namely supernova 1987A. These neutrinos have extensively been used to constrain 
neutrino masses and oscillation parameters @|. 

Neutrinos at much higher energies, above a GeV or so, have not been detected yet due 
to the small fluxes and cross sections. However, the disadvantage of difficult detection is at 
the same time a blessing because it makes these elusive particles reach us unattenuated over 
cosmological distances and from very dense environments where all other particles (except 
gravitational waves) would be absorbed. Such neutrinos are expected to be produced, apart 
from more speculative mechanisms such as decay of superheavy particles, as secondary prod- 
ucts of interactions of the well known cosmic rays in extragalactic and galactic sources as well 
as during propagation. These high energy neutrinos are therefore inextricably linked with 
the physics and astrophysics of cosmic rays. These lectures will therefore discuss high energy 
neutrinos together with cosmic rays. The main focus will thereby be on new possibilities 
opened by new experimental techniques now under consideration or already under construc- 
tion, and resulting probes of new physics at very high energies beyond the electroweak scale. 
We thus begin with a short discussion of the now very active field of ultrahigh energy cosmic 
rays in general. 

After almost 90 years of research on cosmic rays (CRs), their origin is still an open 
question, for which the degree of uncertainty increases with energy: Only below 100 MeV 
kinetic energy, where the solar wind shields protons coming from outside the solar system, 
the sun must give rise to the observed proton flux. The bulk of the CRs up to at least an 
energy of E = 4 x 10 15 eV is believed to originate within our Galaxy. Above that energy, 
which is associated with the so called "knee", the flux of particles per area, time, solid 
angle, and energy, which can be well approximated by broken power laws oc Z? -7 , steepens 
from a power law index 7 ~ 2.7 to one of index ~ 3.2. Above the so called "ankle" at 
E ~ 5 x 10 18 eV, the spectrum flattens again to a power law of index 7 ~ 2.8. This latter 
feature is often interpreted as a cross over from a steeper Galactic component to a harder 
component of extragalactic origin. Fig. |T] shows the measured CR spectrum above 100 MeV, 
up to 3 x 10 20 eV, the highest energy measured so far for an individual CR. 

The conventional scenario assumes that all high energy charged particles are accelerated 
in magnetized astrophysical shocks, whose size and typical magnetic field strength determines 
the maximal achievable energy, similar to the situation in man made particle accelerators. 
The most likely astrophysical accelerators for CR up to the knee, and possibly up to the ankle 
are the shocks associated with remnants of past Galactic supernova explosions, whereas for 
the presumed extragalactic component powerful objects such as active galactic nuclei are 
envisaged. 

The main focus of this contribution will be on ultrahigh energy cosmic rays (UHECRs), 
those with energy > 10 18 eV || 0, ||, [11], [12], see Fig. Q, and neutrinos in the same energy 
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range, which have not been detected yet. For more details on CRs in general the reader is 
referred to recent monographs [111, |ll| . In particular, extremely high energy (EHE)[| cosmic 



rays pose a serious challenge for conventional theories of CR origin based on acceleration 
of charged particles in powerful astrophysical objects. The question of the origin of these 
EHECRs is, therefore, currently a subject of much intense debate and discussions as well 
as experimental efforts; for a short review for the non-experts see Ref. [IB|, for recent brief 



reviews see Ref. ||T7| , |18|j and for a detailed review see Ref. [|19|]. In Sect. 2 we will summarize 



detection techniques and present and future experimental projects. 

The current theories of origin of EHECRs can be broadly categorized into two distinct 
"scenarios": the "bottom-up" acceleration scenario, and the "top-down" decay scenario, 
with various different models within each scenario. As the names suggest, the two scenarios 
are in a sense exact opposite of each other. The bottom- up scenario is just an extension of 
the conventional shock acceleration scenario in which charged particles are accelerated from 
lower energies to the requisite high energies in certain special astrophysical environments. On 
the other hand, in the top-down scenario, the energetic particles arise simply from decay of 
certain sufficiently massive particles originating from physical processes in the early Universe, 
and no acceleration mechanism is needed. 

The problems encountered in trying to explain EHECRs in terms of acceleration mech- 



anisms have been well-documented in a number of studies; see, e.g., Refs. [£(J, [21], |22|1 - Even 
if it is possible, in principle, to accelerate particles to EHECR energies of order 100 EeV in 
some astrophysical sources, it is generally extremely difficult to get the particles come out of 
the dense regions in and/or around the sources without losing much energy. Currently, the 
most favorable sources in this regard are perhaps a class of powerful radio galaxies (see, e.g., 



Refs. |23|, |24| for recent reviews and references to the literature), although the values of the 



relevant parameters required for acceleration to energies > 100 EeV are somewhat on the 



extreme side ||22|| . However, even if the requirements of energetics are met, the main problem 
with radio galaxies as sources of EHECRs is that most of them seem to lie at large cosmo- 
logical distances, 3> 100 Mpc, from Earth. This is a major problem if EHECR particles are 
conventional particles such as nucleons or heavy nuclei. The reason is that nucleons above ~ 
70 EeV lose energy drastically during their propagation from the source to Earth due to the 
Greisen-Zatsepin-Kuzmin (GZK) effect | f25] , 2"B| , namely, photo-production of pions when the 
nucleons collide with photons of the cosmic microwave background (CMB), the mean-free 
path for which is ~ few Mpc [[27]]. This process limits the possible distance of any source 
of EHE nucleons to < 100 Mpc. If the particles were heavy nuclei, they would be photo- 
disintegrated |SR in the CMB and infrared (IR) background within similar distances. 
Thus, nucleons or heavy nuclei originating in distant radio galaxies are unlikely to survive 
with EHECR energies at Earth with any significant flux, even if they were accelerated to 
energies of order 100 EeV at source. In addition, if cosmic magnetic fields are not close to 
existing upper limits, EHECRs are not likely to be deflected strongly by large scale cosmo- 
logical and/or Galactic magnetic fields Thus, EHECR arrival directions should point back to 
their sources in the sky and EHECRs may offer us the unique opportunity of doing charged 



1 We shall use the abbreviation EHE to specifically denote energies E > 10 20 eV, while the abbreviation 
UHE for "Ultra- High Energy" will sometimes be used to denote E > 1 EeV, where 1 EeV = 10 18 eV. Clearly 
UHE includes EHE but not vice versa. 
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particle astronomy. Yet, for the observed EHECR events so far, no powerful sources close 
to the arrival directions of individual events are found within about 100 Mpc p0| , pTf . Very 
recently, it has been suggested by Boldt and Ghosh |3T[ that particles may be accelerated to 
energies ~ 10 21 eV near the event horizons of spinning supermassive black holes associated 
with presently inactive quasar remnants whose numbers within the local cosmological Uni- 
verse (i.e., within a GZK distance of order 50 Mpc) may be sufficient to explain the observed 
EHECR flux. This would solve the problem of absence of suitable currently active sources 
associated with EHECRs. A detailed model incorporating this suggestion, however, remains 
to be worked out. 

There are, of course, ways to avoid the distance restriction imposed by the GZK effect, 
provided the problem of energetics is somehow solved separately and provided one allows 
new physics beyond the Standard Model of particle physics; we shall discuss suggestions 
based on neutrinos in Sect. 3. 

On the other hand, in the top-down scenario, which will be discussed in Sect. 4, the 
problem of energetics is trivially solved from the beginning. Here, the EHECR particles owe 
their origin to decay of some supermassive "X" particles of mass nix ^> 10 20 eV, so that 
their decay products, envisaged as the EHECR particles, can have energies all the way up to 
~ mx- Thus, no acceleration mechanism is needed. The sources of the massive X particles 
could be topological defects such as cosmic strings or magnetic monopoles that could be 
produced in the early Universe during symmetry-breaking phase transitions envisaged in 
Grand Unified Theories (GUTs). In an inflationary early Universe, the relevant topological 
defects could be formed at a phase transition at the end of inflation. Alternatively, the X 
particles could be certain supermassive metastable relic particles of lifetime comparable to or 
larger than the age of the Universe, which could be produced in the early Universe through, 
for example, particle production processes associated with inflation. Absence of nearby 
powerful astrophysical objects such as AGNs or radio galaxies is not a problem in the top- 
down scenario because the X particles or their sources need not necessarily be associated 
with any specific active astrophysical objects. In certain models, the X particles themselves 
or their sources may be clustered in galactic halos, in which case the dominant contribution 
to the EHECRs observed at Earth would come from the X particles clustered within our 
Galactic Halo, for which the GZK restriction on source distance would be of no concern. 

By focusing primarily on "non-conventional" scenarios involving new particle physics be- 
yond the electroweak scale, we do not wish to give the wrong impression that these scenarios 
explain all aspects of EHECRs and UHE neutrinos. In fact, essentially each of the specific 
models that have been studied so far has its own peculiar set of problems. Indeed, the 
main problem of non-astrophysical solutions of the EHECR problem in general is that they 
are highly model dependent. On the other hand, it is precisely because of this reason that 
these scenarios are also attractive — they bring in ideas of new physics beyond the Standard 
Model of particle physics (such as Grand Unification and new interactions beyond the reach 
of terrestrial accelerators) as well as ideas of early Universe cosmology (such as topological 
defects and/or massive particle production in inflation) into the realms of EHECRs where 
these ideas have the potential to be tested by future EHECR experiments. 

The physics and astrophysics of UHECRs is not only intimately linked with the emerging 
field of neutrino astronomy (for reviews see Refs. |3"3]]), but also with the already estab- 



3 



lished field of 7— ray astronomy (for reviews see, e.g., Ref. [j34l) which in turn are important 



sub disciplines of particle astrophysics (for a review see, e.g., Ref. P5[). Indeed, as we shall 
see, all scenarios of UHECR origin, including the top-down models, are severely constrained 
by neutrino and 7— ray observations and limits. In turn, this linkage has important con- 
sequences for theoretical predictions of fluxes of extragalactic neutrinos above a TeV or so 
whose detection is a major goal of next-generation neutrino telescopes (see Sect. 2): If these 
neutrinos are produced as secondaries of protons accelerated in astrophysical sources and 
if these protons are not absorbed in the sources, but rather contribute to the UHECR flux 
observed, then the energy content in the neutrino flux can not be higher than the one in 



UHECRs, leading to the so called Waxman-Bahcall bound [36, 38]. If one of these assump 



tions does not apply, such as for acceleration sources that are opaque to nucleons or in the 
TD scenarios where X particle decays produce much fewer nucleons than 7— rays and neu- 
trinos, the Waxman-Bahcall bound does not apply, but the neutrino flux is still constrained 
by the observed diffuse 7— ray flux in the GeV range. This is true as long as the energy 
fluences produced in 7— rays and neutrinos are comparable, which follows from isospin sym- 
metry if neutrinos are produced by pion production, because 7— rays injected above the pair 
production threshold on the CMB will cascade down to the GeV regime (see Sect. 4.4). 



2 Present and Future UHE CR and Neutrino Experi- 
ments 

The CR primaries are shielded by the Earth's atmosphere and near the ground reveal their 
existence only by indirect effects such as ionization. Indeed, it was the height dependence 
of this latter effect which lead to the discovery of CRs by Hess in 1912. Direct observation 
of CR primaries is only possible from space by flying detectors with balloons or spacecraft. 
Naturally, such detectors are very limited in size and because the differential CR spectrum 
is a steeply falling function of energy (see Fig. [I]), direct observations run out of statistics 
typically around a few 100 TeV. 

Above ~ 100 TeV, the showers of secondary particles created in the interactions of the 
primary CR with the atmosphere are extensive enough to be detectable from the ground. In 
the most traditional technique, charged hadronic particles, as well as electrons and muons 



in these Extensive Air Showers (EAS) are recorded on the ground [^9| with standard instru- 
ments such as water Cherenkov detectors used in the old Volcano Ranch || and Haverah 
Park m experiments, and scintillation detectors which are used now-a-days. Currently op- 
erating ground arrays for UHECR EAS are the Yakutsk experiment in Russia fly]] and the 
Akeno Giant Air Shower Array (AGASA) near Tokyo, Japan, which is the largest one, cov- 
ering an area of roughly 100 km 2 with about 100 detectors mutually separated by about 
lkm [Tj|. The Sydney University Giant Air Shower Recorder (SUGAR) J7| operated until 



1979 and was the largest array in the Southern hemisphere. The ground array technique 
allows one to measure a lateral cross section of the shower profile. The energy of the shower- 
initiating primary particle is estimated by appropriately parametrizing it in terms of a mea- 
surable parameter; traditionally this parameter is taken to be the particle density at 600 m 
from the shower core, which is found to be quite insensitive to the primary composition and 
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the interaction model used to simulate air showers. 

The detection of secondary photons from EAS represents a complementary technique. 
The experimentally most important light sources are the fluorescence of air nitrogen excited 
by the charged particles in the EAS and the Cherenkov radiation from the charged particles 
that travel faster than the speed of light in the atmospheric medium. The first source is 
practically isotropic whereas the second one produces light strongly concentrated on the 
surface of a cone around the propagation direction of the charged source. The fluorescence 
technique can be used equally well for both charged and neutral primaries and was first used 
by the Fly's Eye detector [12| and will be part of several future projects on UHECRs (see 
below). The primary energy can be estimated from the total fluorescence yield. Information 
on the primary composition is contained in the column depth X max (measured in gem -2 ) at 
which the shower reaches maximal particle density. The average of X max is related to the 
primary energy E by 

(X max > = X^ In (^L) . (1) 



Here, X' Q is called the elongation rate and Eq is a characteristic energy that depends on 
the primary composition. Therefore, if X max and X' Q are determined from the longitudinal 
shower profile measured by the fluorescence detector, then Eq and thus the composition, can 
be extracted after determining the energy E from the total fluorescence yield. Comparison 
of CR spectra measured with the ground array and the fluorescence technique indicate 
systematic errors in energy calibration that are generally smaller than ~ 40%. For a more 
detailed discussion of experimental EAS analysis with the ground array and the fluorescence 
technique see, e.g., Refs. |f40f . 

As an upscaled version of the old Fly's Eye Cosmic Ray experiment, the High Resolution 
Fly's Eye detector has started to take data at Utah, USA |42| . Taking into account a duty 
cycle of about 10% (a fluorescence detector requires clear, moonless nights), the effective 
aperture of this instrument will be ~ 350(1000) km 2 sr at 10(100) EeV, on average about 6 
times the Fly's Eye aperture, with a threshold around 10 17 eV. Another project utilizing the 
fluorescence technique is the Japanese Telescope Array (^3j which is currently in the proposal 
stage. If approved, its effective aperture will be about 10 times that of Fly's Eye above 
10 17 eV, and it would also be used as a Cherenkov detector for TeV 7— ray astrophysics. The 
largest project presently under construction is the Pierre Auger Giant Array Observatory [ pH 



planned for two sites, one in Argentina and another in the USA for maximal sky coverage. 
Each site will have a 3000 km ground array. The southern site will have about 1600 particle 
detectors (separated by 1.5 km each) overlooked by four fluorescence detectors. The ground 
arrays will have a duty cycle of nearly 100%, leading to an effective aperture about 30 times 
as large as the AGASA array. The corresponding cosmic ray event rate above 10 20 eV will 
be about 50 events per year. About 10% of the events will be detected by both the ground 
array and the fluorescence component and can be used for cross calibration and detailed EAS 
studies. The energy threshold will be around 10 18 eV, with full sensitivity above 10 19 eV. 

Recently NASA initiated a concept study for detecting EAS from space [[4l| f|7| by ob- 
serving their fluorescence light from an Orbiting Wide-angle Light-collector (OWL). This 
would provide an increase by another factor ~ 50 in aperture compared to the Pierre Auger 
Project, corresponding to an event rate of up to a few thousand events per year above 10 20 eV. 
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Similar concepts such as the Extreme Universe Space Observatory (EUSO) which is part 
of the AirWatch program [|9| and of which a prototype may be tested on the International 
Space Station are also being discussed. It is possible that the OWL and AirWatch efforts will 
merge. The energy threshold of such instruments would be between 10 19 and 10 20 eV. This 
technique would be especially suitable for detection of very small event rates such as those 
caused by UHE neutrinos which would produce deeply penetrating EAS (see Sect. 3.2). For 
more details on these recent experimental considerations see Ref. 
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High energy neutrino astronomy is aiming towards a kilometer scale neutrino observatory. 
The major technique is the optical detection of Cherenkov light emitted by muons created 
in charged current reactions of neutrinos with nucleons either in water or in ice. The largest 
pilot experiments representing these two detector media are the now defunct Deep Undersea 
Muon and Neutrino Detection (DUMAND) experiment |5l| in the deep sea near Hawai and 
the Antarctic Muon And Neutrino Detector Array (AMANDA) experiment |52[] in the South 
Pole ice. Another water based experiment is situated at Lake Baikal ||53|| . Next generation 
deep sea projects include the French Astronomy with a Neutrino Telescope and Abyss envi- 
ronmental RESearch (ANTARES) |54[] and the underwater Neutrino Experiment SouthwesT 
Of GReece (NESTOR) project in the Mediterranean [Q, whereas ICECUBE represents 
the planned kilometer scale version of the AMANDA detector. Also under consideration are 
neutrino detectors utilizing techniques to detect the radio pulse from the electromagnetic 
showers created by neutrino interactions in ice and other materials (see Ref. |57|]). This 
technique could possibly be scaled up to an effective area of 10 4 km 2 and a prototype is rep- 



resented by the Radio Ice Cherenkov Experiment (RICE) experiment at the South Pole [58 
The Goldstone radio telescope has already put an upper limit on UHE neutrino fluxes 



from the non-observation of radio pulses from showers induced by neutrinos interacting in 
the moons rim The sensitivity of existing underwater acoustic arrays to neutrino in- 
duced showers is also being studied [^]. Furthermore, neutrinos can also initiate horizontal 
EAS which can be detected by giant ground arrays and fluorescence detectors such as the 
Pierre Auger Project U63 |. It has also been shown that the sensitivity of such experiments to 
UHE neutrinos could be significantly enhanced by triggering on Earth-skimming showers, i.e. 
showers with zenith angles slightly larger than 90° |>|, |65], |66|, especially for r— neutrinos. 



Finally, as mentioned above, deeply penetrating EAS could be detected from space by in- 
struments such as the proposed space based AirWatch type detectors |E6l [|7], 





49|. 


. 32 


,67 



and 



some recent overviews on neutrino astronomy can be found in Ref. [33 



3 Neutrino Interactions and Propagation 
3.1 Neutrino propagation 

The propagation of UHE neutrinos is governed mainly by their interaction with the relic 
neutrino background (RNB). In this section we give a short overview over the relevant 
interactions within the general framework for particle propagation used in the present con- 
tribution. 

The average squared CM energy for interaction of an UHE neutrino of energy E with a 
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relic neutrino of energy e is given by 



(s) ~ (45 GeV) 1 



E 



10- 3 eV/ V10 15 GeV 



(2) 



If the relic neutrino is relativistic, then e ~ 3^(1+^/4) in Eq. (0), where T v ~ 1.9(l+z) K = 
1.6 x 10~ 4 (1 + z) eV is the temperature at redshift z and r/b < 50 is the dimensionless 
chemical potential of relativistic relic neutrinos. For nonrelativistic relic neutrinos of mass 



m v < 20 eV, e ~ max [3T U 



rriv 



Note that Eq. (0) implies interaction energies that are 



typically smaller than electroweak energies even for UHE neutrinos, except for energies near 
the Grand Unification scale, E > 10 15 GeV, or if m u > 1 eV. In this energy range, the 
cross sections are given by the Standard Model of electroweak interactions which are well 
confirmed experimentally. Physics beyond the Standard Model is, therefore, not expected to 
play a significant role in UHE neutrino interactions with the low energy relic backgrounds. 

The dominant interaction mode of UHE neutrinos with the RNB is the exchange of a 
W ± boson in the t-channel (z/j + Vj — > Zj + lj), or of a Z° boson in either the s-channel 
(z/j + 9i — > //) or the t-channel (z^ + Vj — > i/j + Vj) |68], [69|, |70], |71fl . Here, i, j stands for either 
the electron, muon, or tau flavor, where i / j for the first reaction, I denotes a charged 
lepton, and / any charged fermion. If the latter is a quark, it will, of course, subsequently 
fragment into hadrons. As an example, the differential cross section for s-channel production 
of Z° is given by 



da 



dfi 



GpS 
An 



Ml 



[s - miy + Mfr| 



0£(i + aO 2 + s£(i-aO 5 



(3) 



where Gf is the Fermi constant, Mz and V z are mass and lifetime of the Z°, gi and gn are 
the usual dimensionless left- and right-handed coupling constants for /, and \i* is the cosine 
of the scattering angle in the CM system. 

The t-channel processes have cross sections that rise linearly with s up to s 
M w the W ± mass, above which they are roughly constant with a value a t (s > M w ) 
G F M^r ~ 10~ 34 cm 2 . Using Eq. (§) this yields the rough estimate 



Myy, With 



a t (E,e) 



~ mm 



~ mm 



10" 34 ,10 



-li 



MeV z 
10" 34 ,3 x 10~ 39 



cm 



(4) 



( E 



10" 3 eV/ V10 20 eV 



cm 2 . 



In contrast, within the Standard Model the neutrino-nucleon cross section roughly behaves 
as 

a uN (E) ~ 10- 31 ( J E/10 20 eV) ' 4 cm 2 (5) 

for E > 10 15 eV (see Eq.(^) below). Interactions of UHE neutrinos with nucleons are, how- 
ever, still negligible compared to interactions with the RNB because the RNB particle density 
is about ten orders of magnitude larger than the baryon density. The only exception could 
occur near Grand Unification scale energies and at high redshifts and/or if contributions 
to the neutrino-nucleon cross section from physics beyond the Standard Model dominate at 
these energies (see Sect. 3.3 below). 
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It has recently been pointed out JF^I that above the threshold for W ± production the 
process v + 7 — > IW + becomes comparable to the vv processes discussed above. Fig. |3] 
compares the cross sections relevant for neutrino propagation at CM energies around the 
electroweak scale. Again, for UHE neutrino interactions with the RNB the relevant CM 
energies can only be reached if (a) the UHE neutrino energy is close to the Grand Unification 
scale, or (b) the RNB neutrinos have masses in the eV regime, or (c) at redshifts z > 10 3 . 
Even then the z/7 process never dominates over the vv process. 

At lower energies there is an additional z/7 interaction that was recently discussed as 

potentially important besides the vv processes: Using an effective Lagrangian derived from 

the Standard Model, Ref. J7| obtained the result a 7+ „^ 7+7+ „(s) ~ 9 x 1(T 56 (s/MeV 2 ) 5 cm 2 , 

2 



supposed to be valid at least up to s < 10 MeV . Above the electron pair production 
threshold the cross section has not been calculated because of its complexity but is likely 
to level off and eventually decrease. Nevertheless, if the s 5 behavior holds up to s ~ a few 
hundred MeV 2 , comparison with Eq. (|J) shows that the process 7 + z/ — > 7 + 7 + z/ would 
start to dominate and influence neutrino propagation around E ~ 3 x 10 17 (e/10 -3 eV) eV, 



as was pointed out in Ref. 74 



For a given source distribution, the contribution of the "direct" neutrinos to the flux can 
be computed from (t is the age of the Universe) 

j(E) ~ —t / d Zi (l + ^)- 9/2 $[(l + Zi)E, Zi ] , (6) 

o7T JO 

up to the interaction redshift z(E), i.e. the average redshift from which a neutrino of present 
day energy E could have propagated without interacting. This approximation neglects the 
secondary neutrinos and the decay products of the leptons created in the neutral current 
and charged current reactions of UHE neutrinos with the RNB discussed above. Similarly 
to the EM case, these secondary particles can lead to neutrino cascades developing over 
cosmological redshifts |7D[ . 

Approximate expressions for the interaction redshift for the processes discussed above 



have been given in Refs. p4| , [F5[ for CM energies below the electroweak scale, assuming rel- 
ativists, nondegenerate relic neutrinos, m u < T u , and rjb <C 1. Approaching the electroweak 
scale, a resonance occurs in the interaction cross section for s-channel Z° exchange at the Z° 
mass, s = Mf ~ (91 GeV) 2 , see Eq. ([IJ). The absorption redshift for the corresponding neu- 
trino energy, E ~ 10 15 GeV(e/10 -3 eV) -1 drops to a few (or less for a degenerate, relativistic 
RNB) and asymptotically approaches constant values of a few tens at higher energies. 

An interesting situation arises if the RNB consists of massive neutrinos with m v ~ 
1 eV: Such neutrinos would constitute hot dark matter which is expected to cluster [f76| , 
for example, in galaxy clusters. This would potentially increase the interaction probability 
for any neutrino of energy within the width of the Z° resonance at E = M^/2m u = 4 x 
10 21 (eV/m l/ ) eV. Recently it has been suggested that the stable end products of the "Z- 
bursts" that would thus be induced at close-by distances (< 50 Mpc) from Earth may explain 
the highest energy cosmic rays |T7], and may also provide indirect evidence for neutrino 
hot dark matter. These end products would be mostly nucleons and 7— rays with average 
energies a factor of ~ 5 and ~ 40 lower, respectively, than the original UHE neutrino. As a 
consequence, if the UHE neutrino was produced as a secondary of an accelerated proton, the 
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energy of the latter would have to be at least a few 10 22 eV |77| , making Z-bursts above GZK 
energies more likely to play a role in the context of non-acceleration scenarios. Moreover, it 
has subsequently been pointed out |79fl that Z production is dominated by annihilation on 



the non-clustered massive RNB compared to annihilation with neutrinos clustering in the 
Galactic halo or in nearby galaxy clusters. As a consequence, for a significant contribution of 
neutrino annihilation to the observed EHECR flux, a new class of neutrino sources, unrelated 
to UHECR sources, seems necessary. This has been confirmed by more detailed numerical 
simulations 



80 



where it has, however, also been demonstrated that the most significant 
contribution could come from annihilation on neutrino dark matter clustering in the Local 
Supercluster by amounts consistent with expectations. In the absence of any assumptions 
on the neutrino sources, the minimal constraint comes from the unavoidable production of 
secondary 7— rays contributing to the diffuse flux around 10 GeV measured by EGRET: If 
the Z-burst decay products are to explain EHECR, the massive neutrino overdensity f v over 
a length scale l u has to satisfy f u > 20 {l v /h Mpc)^ 1 , provided that only neutrinos leave the 
source, a situation that may arise in top-down models if the X particles decay exclusively 
into neutrinos (see Fig. |8| below for a model involving topological defects and Ref. for a 
scenario involving decaying superheavy relic particles). If, instead, the total photon source 
luminosity is comparable to the total neutrino luminosity, as in most models, the EGRET 
constraint translates into the more stringent requirement f u > 10 3 (/^/5 Mpc) -1 . This bound 
can only be relaxed if most of the EM energy is radiated in the TeV range where the Universe 
is more transparent fl80 |. Fig. |] shows an example of this situation. 



Furthermore, the Z-burst scenario requires sources that are optically thick for acceler- 
ated protons with respect to photo-pion production because otherwise the observable proton 
flux below the GZK cutoff would be comparable to the neutrino flux p9|, |79|. A systematic 



parameter study of the required overdensity, based on analytical flux estimates, has been 
performed in Ref. [90|. Recently it has been noted that a degenerate relic neutrino back- 



ground would increase the interaction probability and thereby make the Z-burst scenario 

A neutrino asymmetry of order unity is not excluded phenomenolog- 



91 



more promising 

ically [[52|] and can be created in the early Universe, for example, through the Affleck-Dine 
baryogenesis mechanism or due to neutrino oscillations. The authors of Ref. |9l[] pointed 
out that for a neutrino mass ra v ~ 0.07 eV, a value suggested by the Super-Kamiokande ex- 
periment || , and for sources at redshifts of a few, the flux of secondary Z-decay products is 
maximal for a RNB density parameter Q u ~ 0.01. Such neutrino masses, however, require 
the sources to produce neutrinos at least up to 10 22 eV. 

UHE neutrinos from the decay of pions, that are produced by interactions of accelerated 
protons in astrophysical sources, must have originated within redshifts of a few. Moreover, in 
most conventional models their flux is expected to fall off rapidly above 10 20 eV. Examples 
are production in active galactic nuclei within hadronic models |95], [97|, |98], p9| , 
and "cosmogenic" neutrinos from interactions of UHECR nucleons 
cutoff) with the CMB (see, e.g.,Refs. |TDg, pIH 



near or above the GZK 
The latter source is the only one that is 
guaranteed to exist due to existence of UHECRs near the GZK cutoff, but the fluxes are 
generally quite small. Therefore, interaction of these UHE neutrinos with the RNB, that 
could reveal the latter's existence, can, if at all, be important only if the relic neutrinos have 
a mass m u > 1 eV . Due to the continuous release of UHE neutrinos up to much higher 
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redshifts, most top-down scenarios would imply substantially higher fluxes that also extend 
to much higher energies |75|] (see Sect. 4.4 below). Certain features in the UHE neutrino 
spectrum predicted within such top-down scenarios, such as a change of slope for massless 
neutrinos [|70| or a dip structure for relic neutrino masses of order 1 e V |7l], [78[] , have therefore 
been proposed as possibly the only way to detect the RNB. However, some of the scenarios 
at the high end of neutrino flux predictions have recently been ruled out based on constraints 
on the accompanying energy release into the EM channel. 

The recent claim that the Z-burst mechanism can explain the EHECR flux even without 
overdensity, f u = 1, and even allows to determine the neutrino mass m v ~ leV ||102|| was 
based on assumptions that we characterized here as highly unrealistic: Sources accelerating 
nuclei to > 10 23 eV while being completely opaque to both the primaries and the secondary 
photons. The energy fluence of the latter, however, after possibly being reprocessed to lower 
energies by EM cascading, must be comparable to the neutrino energy fluence by simple 
isospin symmetry in the production of charged and neutral pions. 

Since in virtually all models UHE neutrinos are created as secondaries from pion de- 
cay, i.e. as electron or muon neutrinos, r— neutrinos can only be produced by a flavor 
changing W t-channel interaction with the RNB. The flux of UHE r— neutrinos is there- 
fore usually expected to be substantially smaller than the one of electron and muon neu- 
trinos, if no neutrino oscillations take place at these energies. However, the recent ev- 
idence from the Superkamiokande experiment for nearly maximal mixing between muon 
and t— neutrinos with |Am 2 | = |m 2 — m 2 J ~ 5 x 10~ 3 eV 2 would imply an oscil- 
lation length of L osc = 2E/\Am 2 \ = 2.6 x KT 6 (£/PeV)(| Am 2 |/5 x 10" 3 eV 2 )- 1 pc and, 
therefore, a rough equilibration between muon and r— neutrino fluxes from any source at 
a distance larger than L osc ||104|| . Turning this around, one sees that a source at distance 



d emitting neutrinos of energy E is sensitive to neutrino mixing with | Am 2 1 = 2E/d ~ 
1.3 x 1CU 16 (E/PeV)(d/100 Mpc)" 1 eV 2 ||105 , 106 . Under certain circumstances, resonant 



conversion in the potential provided by the RNB clustering in galactic halos may also influ- 
ence the flavor composition of UHE neutrinos from extraterrestrial sources ||107 |. In addition, 



such huge cosmo logical baselines can be sensitive probes of neutrino decay |108 [ 



3.2 Neutrino detection 

We now turn to a discussion of UHE neutrino interactions with matter relevant for neutrino 
detection. UHE neutrinos can be detected by detecting the muons produced in ordinary 
matter via charged-current reactions with nucleons; see, e.g., Refs. [|109| , |110| , |111|| for recent 



discussions. Corresponding cross sections are calculated by folding the fundamental stan- 
dard model quark-neutrino cross section with the distribution function of the partons in the 
nucleon. These cross sections are most sensitive to the abundance of partons of fractional 
momentum x ~ My V /2rriNE, where E is the neutrino energy. For the relevant squared mo- 
mentum transfer, Q 2 ~ M^, these parton distribution functions have been measured down 
to x ~ 0.02 [ 112|| . (It has been suggested that observation of the atmospheric neutrino flux 



with future neutrino telescopes may probe parton distribution functions at much smaller x 
currently inaccessible to colliders ||113|| ). Currently, therefore, neutrino-nucleon cross sections 



for E > 10 14 eV can be obtained only by extrapolating the parton distribution functions to 
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lower x. Above 10 19 eV, the resulting uncertainty has been estimated to be a factor 2 ||1 10 



whereas within the dynamical radiative parton model it has been claimed to be at most 
20 % An intermediate estimate using the CTEQ4-DIS distributions can roughly be 



parameterized by |[L10|| 

a uN (E) ~ 2.36 x 10~ 32 (£/10 19 eV) a363 cm 2 (10 16 eV < E < 10 21 eV) . (7) 
Improved calculations including non- leading logarithmic contributions in 1/x have recently 



been performed in Ref. 114 . The results for the neutrino-nucleon cross section differ by less 



than a factor 1.5 with Refs. ||1 10| , |111|| even at 10 21 eV. 



However, more recently it has been argued that the neutrino-nucleon cross-section calcu- 
lated within the Standard Model becomes unreliable for E > 2 x 10 17 eV: the authors of ||1 1 5| 



used the 0(g 2 ) expression for the elastic forward scattering amplitude to derive via the optical 
theorem the bound <j v n < 9.3 x 10 _33 cm 2 for the total cross-section. Current parton distri- 
bution functions (pdf's) predict a violation of this unitarity bound above E > 2 x 10 17 eV. 
The authors of [|115|| argue that the large 0(g A ) corrections to the forward scattering ampli- 
tude necessary to restore unitarity signal a breakdown of electroweak perturbation theory. 
Alternatively, large changes in the evolution of the parton distribution functions have to set 
in soon after the kinematical range probed by HERA. A large 0(g A ) correction to the for- 
ward scattering amplitude is, however, not surprising because the 0(g 2 ) amplitude contains 
no resonant contribution and is real. In particular, the total cross-section is therefore not 
only bounded by a constant but zero at 0(g 2 ), and the imaginary part of the box diagram 
of 0(g 4 ) is the first contribution to the total neutrino-nucleon cross-section ||116|| . 

Interestingly, it has been shown that the increasing target mass provided by the Earth for 
increasing zenith angles below the horizontal implies that the rate of up-going air showers in 
UHECR detectors does not decrease with decreasing neutrino-nucleon cross section but may 
even increase f6~6fl . Thus, cross sections smaller than Eq. (0) do not lead to reduced event 



rates in UHECR detectors and can be measured from the angular distribution of events. 
UHECR and neutrino experiments can thus contribute to measure cross sections at energies 
inaccessible in accelerator experiments! 

Neutral-current neutrino-nucleon cross sections are expected to be a factor 2-3 smaller 
than charged-current cross sections at UHE and interactions with electrons only play a role 
at the Glashow resonance, v e e — > W, at E = 6.3 x 10 15 eV. Furthermore, cross sections of 
neutrinos and anti-neutrinos are basically identical at UHE. Radiative corrections influence 
the total cross section negligibly compared to the parton distribution uncertainties, but may 
lead to an increase of the average inelasticity in the outgoing lepton from ~ 0.19 to ~ 0.24 
at E ~ 10 20 eV |117|, although this would probably hardly influence the shower character. 



Neutrinos propagating through the Earth start to be attenuated above ~ 100 TeV due to 
the increasing Standard Model cross section as indicated by Eq. (^). Detailed integrations 
of the relevant transport equations for muon neutrinos above a TeV have been presented in 



Ref. ||114|| , and, for a general cold medium, in Ref. ||118|| . In contrast, r— neutrinos with en- 
ergy up to ~ 100 PeV can penetrate the Earth due to their regeneration from r decays [ 106 . 
As a result, a primary UHE r— neutrino beam propagating through the Earth would cas- 
cade down below ~ 100 TeV and in a neutrino telescope could give rise to a higher total 
rate of upgoing events as compared to downgoing events for the same beam arriving from 
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above the horizon. As mentioned above, a primary r— neutrino beam could arise even in 
scenarios based on pion decay, if — v T mixing occurs with the parameters suggested by the 
Super-Kamiokande results ||104| . In the PeV range, r— neutrinos can produce characteristec 



"double-bang" events where the first bang would be due to the charged-current production 
by the r— neutrino of a r whose decay at a typical distance ~ 100 m would produce the 
second bang |105|| . These effects have also been suggested as an independent astrophysical 



test of the neutrino oscillation hypothesis. In addition, isotropic neutrino fluxes in the en- 
ergy range between 10 TeV and 10 PeV have been suggested as probes of the Earth's density 
profile, whereby neutrino telescopes could be used for neutrino absorption tomography ||119| . 



3.3 New Interactions 

It has been suggested that the neutrino-nucleon cross section, <j u n, can be enhanced by new 
physics beyond the electroweak scale in the center of mass (CM) frame, or above about a 
PeV in the nucleon rest frame ||120|, |121|, |122j|. Neutrino induced air showers may therefore 



rather directly probe new physics beyond the electroweak scale. 

The lowest partial wave contribution to the cross section of a point-like particle is con- 



strained by unitarity to be not much larger than a typical electroweak cross section ||123 
However, at least two major possibilities allowing considerably larger cross sections have 
been discussed in the literature for which unitarity bounds need not be violated. In the first, 
a broken SU(3) gauge symmetry dual to the unbroken SU(3) color gauge group of strong 
interaction is introduced as the "generation symmetry" such that the three generations of 
leptons and quarks represent the quantum numbers of this generation symmetry. In this 
scheme, neutrinos can have close to strong interaction cross sections with quarks. In ad- 
dition, neutrinos can interact coherently with all partons in the nucleon, resulting in an 
effective cross section comparable to the geometrical nucleon cross section. This model lends 



itself to experimental verification through shower development altitude statistics ||120|| . 

The second possibility consists of a large increase in the number of degrees of freedom 
above the electroweak scale ||124|| . A specific implementation of this idea is given in theories 



with n additional large compact dimensions and a quantum gravity scale M 4+n ~ TeV that 
has recently received much attention in the literature [ |125|| because it provides an alternative 



solution (i.e., without supersymmetry) to the hierarchy problem in grand unifications of 
gauge interactions. The cross sections within such scenarios have not been calculated from 
first principles yet. Within the field theory approximation which should hold for squared 



CM energies s < Mf +n , the spin 2 character of the graviton predicts a g ~ s 2 /M| +n ||126|| . For 
s ^> M| +n , several arguments based on unitarity within field theory have been put forward. 
The emission of massive Kaluza-Klein (KK) graviton modes associated with the increased 
phase space due to the extra dimensions leads to the rough estimate (derived for n = 2) ||126 



a n ~ — ~ 10" 27 — - - cm 2 



9 ~ Mt +n ~ \ TeV J Vl0 20 eW 



where in the last expression we specified to a neutrino of energy E hitting a nucleon at rest. 
A more detailed calculation taking into account scattering on individual partons leads to 



similar orders of magnitude ||122 |. Note that a neutrino would typically start to interact in 
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the atmosphere for a uN > 10 



-27 cm 2 



i.e. in the case of Eq. (|[) for E > 10 20 eV, assuming 



M 4+n ~ 1 TeV. For cross sections such large the neutrino therefore becomes a primary 
candidate for the observed EHECR events. However, since in a neutral current interaction 
the neutrino transfers only about 10% of its energy to the shower, the cross section probably 
has to be at least a few 10~ 26 cm 2 to be consistent with observed showers which start within 
the first 50 gem -2 of the atmosphere [ |127| , |128| ]. A specific signature of this scenario would 
be the absence of any events above the energy where o g grows beyond ~ 10~ 27 cm 2 in 
neutrino telescopes based on ice or water as detector medium ||129| , and a hardening of the 



spectrum above this energy in atmospheric detectors such as the Pierre Auger Project |44 



and the proposed space based AirWatch type detectors [[46], |48|, [49 
to Eq. i 



Furthermore, according 



I, the average atmospheric column depth of the first interaction point of neutrino 
induced EAS in this scenario is predicted to depend linearly on energy. This should be 
easy to distinguish from the logarithmic scaling expected for nucleons, nuclei, and 7— rays. 
To test such scalings one can, for example, take advantage of the fact that the atmosphere 
provides a detector medium whose column depth increases from ~ 1000 g/cm 2 towards the 
zenith to ~ 36000 g/cm 2 towards horizontal arrival directions. This probes cross sections 
in the range ~ 10 -29 — 10~ 27 cm 2 . Due to the increased column depth, water/ice detectors 



10~ 31 - 10 



-29 cm 2 



301 which could be relevant 



would probe cross sections in the range 
for TeV scale gravity models [ |131|| . 

From a perturbative point of view within string theory, a U N can be estimated as follows: 
Individual amplitudes are expected to be suppressed exponentially above the string scale 
M s which again for simplicity we assume here to be comparable to M4 +n . This can be 
interpreted as a result of the finite spatial extension of the string states. In this case, the 
neutrino nucleon cross section would be dominated by interactions with the partons carrying 
a momentum fraction x ~ M 2 /s, leading to g27| 



o vN ~ ^ln( S /M 2 )( S /M s 2 )°- 363 ^6xl0- 29 (^ 



-4.726 



E 



10 20 eV 



u.:!(i:; 



X 



1 + 0.08 In 



E 



10 20 eV 



-0.16 In (^L) 



cm 



(9) 



This is probably too small to make neutrinos primary candidates for the highest energy show- 
ers observed, given the fact that complementary constraints from accelerator experiments 
result in M s > 1 TeV ||132 ]. On the other hand, general arguments on the production of 
"string balls" or small black holes from two point particles represented by light strings ||133|| 
leads to the asymptotic scaling a ~ (s/M 2 ) 1 /*™ 4-1 ) for s > M 2 for the fundamental neutrino- 
parton cross section. This could lead to values for a U N larger than Eq. ||134j| . Some other 
recent work seems to imply that in TeV string models cross sections, if not sufficient to make 
neutrinos UHECR primary candidates, could at least be significantly larger than Standard 
Model cross sections ||135||. Thus, an experimental detection of the signatures discussed in 



this section could lead to constraints on some string-inspired models of extra dimensions. 

There are, however, severe astrophysical and cosmological constraints on M 4+ „ which 
result from limiting the emission of bulk gravitons into the extra dimensions. The strongest 
constraints in this regard come from the production due to nucleon-nucleon bremsstrahlung 



in type II supernovae [136] and their subsequent decay into a diffuse background of 7— rays in 
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the MeV range fll37fl . The latter read M 6 > 84TeV, M 7 > 7TeV, for n = 2,3, respectively, 
and, therefore, n > 5 is required if neutrino primaries are to serve as a primary candidate for 
the EHECR events observed above 10 20 eV. This assumes a toroidal geometry of the extra 
dimensions with equal radii given by 

where Mp\ denotes the Planck mass. The above lower bounds on M^ +n thus translate into 
the corresponding upper bounds r 2 < 0.9 x lCr 4 mm, r 3 < 0.19 x 10~ 6 mm, respectively. 
Still stronger but somewhat more model dependent bounds result from the production of 
KK modes during the reheating phase after inflation. For example, a bound M$ > 500 TeV 
has been reported ||138|| based on the contribution to the diffuse 7— ray background in the 100 



MeV region. We note, however, that all astrophysical and cosmological bounds are changed 
in more complicated geometries of extra dimensions [ |189[ . 



The neutrino primary hypothesis of EHECR together with other astrophysical and cos- 
mological constraints thus provides an interesting testing ground for theories involving large 
compact extra dimensions representing one possible kind of physics beyond the Standard 
Model. In this context, we mention that in theories with large compact extra dimensions 
mentioned above, Newton's law of gravity is expected to be modified at distances smaller 



than the length scale given by Eq. (10). Indeed, there are laboratory experiments measur- 
ing gravitational interaction at small distances (for a recent review of such experiments see 
Ref. [|140]| ) , which also probe these theories. Thus, future EHECR experiments and gravita- 
tional experiments in the laboratory together have the potential of providing rather strong 
tests of these theories. These tests would be complementary to constraints from collider 
experiments |[L32|| . 

Independent of theoretical arguments, the EHECR data can be used to put constraints 
on cross sections satisfying <j u n(E > 10 19 eV) < 10~ 27 cm 2 . Particles with such cross sections 
would give rise to horizontal air showers. The Fly's Eye experiment established an upper 



limit on horizontal air showers 38 . The non-observation of the neutrino flux expected from 



pions produced by EHECRs interacting with the CMB the results in the limit [fL41| , |130[ | 



a^(10 17 eV) < lxKT 29 /?/ 1 / 2 
<7„tf(10 18 eV) < SxKT 30 /^/ 2 



cm 



2 



cm 



2 



a^(10 19 eV) < 5x lO-^/y^cm 2 , (11) 

where y is the average energy fraction of the neutrino deposited into the shower (y = 1 
for charged current reactions and y ~ 0.1 for neutral current reactions). Neutrino fluxes 
predicted in various scenarios are shown in Fig. ||]. The projected sensitivity of future ex- 
periments such as the Pierre Auger Observatories and the AirWatch type satellite projects 
indicate that the cross section limits Eq. (|TTD could be improved by up to four orders of mag- 
nitude, corresponding to one order of magnitude in M s or M 4+n . This would close the window 
between cross sections allowing horizontal air showers, <7 v n{E 10 19 eV) < 10 _27 cm 2 , and 
the Standard Model value Eq. (|?p. 
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We note in this context that, only assuming 3 + 1 dimensional field theory, consistency of 
the UHE vN cross section with data at electroweak energies does not lead to very stringent 
constraints: Relating the cross section to the uN elastic amplitude in a model independent 
way yields ||145|| 



a(E) < 3 x 10 



-24 



cm 



(12) 



33 cm 2 



,10 19 eV 

However, it has been argued | |115| | (see discussion above), cross sections <j v n > 9.3x10 
could signal a breakdown of perturbation theory. 

In the context of conventional astrophysical sources, the relevant UHE neutrino primaries 
could, of course, only be produced as secondaries in interactions with matter or with low 
energy photons of protons or nuclei accelerated to energies of at least 10 21 eV. This implies 
strong requirements on the possible sources. 

In addition, neutrino primaries with new interactions would predict a significant correla- 
tion of UHECR arrival directions with high redshift objects. Indeed, various claims recently 
occured in the literature for significant angular correlations of UHECRs with certain astro- 
physical objects at distances too large for the primaries to be nucleons, nuclei, or 7— rays. 
Farrar and Biermann reported a possible correlation between the arrival direction of the five 



highest energy CR events and compact radio quasars at redshifts between 0.3 and 2.2 |146 



Undoubtedly, with the present amount of data the interpretation of such evidence for a cor- 
relation remains somewhat subjective, as is demonstrated by the criticism of the statistical 
analysis in Ref. | 146 | by Hoffman | 147| | and the reply by Farrar and Biermann | 148 |). Also, a 



new analysis with the somewhat larger data set now available did not support such correla- 
tions ||149|| . This is currently disputed since another group claims to have found a correlation 
on the 99.9% confidence level [|150j| . Most recently, a correlation between UHECRs of energy 
E > 4 x 10 19 eV and BL Lacertae objects at redshifts z > 0.1 was claimed | |151| |. None of 
these claims are convincing yet but confirmation or refutation should be possible within the 
next few years by the new experiments. Clearly, a confirmation of one of these correlations 
would be exciting as it would probably imply new physics such as neutrinos with new inter- 
actions, as yet unknown, e.g. supersymmetric, particles, or even more radical propositions 
such as violation of Lorentz invariance (for a discussion of these latter possibilities see, e.g., 
Refs. HI 



4 Top-Down Scenarios 
4.1 The Main Idea 

As mentioned in the introduction, all top-down scenarios involve the decay of X particles of 
mass close to the GUT scale which can basically be produced in two ways: If they are very 
short lived, as usually expected in many GUTs, they have to be produced continuously. The 
only way this can be achieved is by emission from topological defects left over from cosmo- 
logical phase transitions that may have occurred in the early Universe at temperatures close 
to the GUT scale, possibly during reheating after inflation. Topological defects necessarily 
occur between regions that are causally disconnected, such that the orientation of the or- 
der parameter associated with the phase transition can not be communicated between these 
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regions and consequently will adopt different values. Examples are cosmic strings (similar 
to vortices in superffuid helium), magnetic monopoles, and domain walls (similar to Bloch 
walls separating regions of different magnetization in a ferromagnet). The defect density is 
consequently given by the particle horizon in the early Universe and their formation can even 
be studied in solid state experiments where the expansion rate of the Universe corresponds 
to the quenching speed with which the phase transition is induced [152|. The defects are 



topologically stable, but in the cosmological case time dependent motion leads to the emis- 
sion of particles with a mass comparable to the temperature at which the phase transition 
took place. The associated phase transition can also occur during reheating after inflation. 

Alternatively, instead of being released from topological defects, X particles may have 
been produced directly in the early Universe and, due to some unknown symmetries, have a 
very long lifetime comparable to the age of the Universe. In contrast to Weakly-Interacting 
Massive Particles (WIMPS) below a few hundred TeV which are the usual dark matter can- 
didates motivated by, for example, supersymmetry and can be produced by thermal freeze 
out, such superheavy X particles have to be produced non-thermally. Several such mech- 
anisms operating in the post- inflationary epoch in the early Universe have been studied. 
They include gravitational production through the effect of the expansion of the background 
metric on the vacuum quantum fluctuations of the X particle field, or creation during re- 
heating at the end of inflation if the X particle field couples to the inflaton field. The latter 
case can be divided into three subcases, namely "incoherent" production with an abundance 
proportional to the X particle annihilation cross section, non-adiabatic production in broad 
parametric resonances with the oscillating inflaton field during preheating (analogous to 
energy transfer in a system of coupled pendula), and creation in bubble wall collisions if 
inflation is completed by a first order phase transition. In all these cases, such particles, 
also called 'WIMPZILLAs", would contribute to the dark matter and their decays could 
still contribute to UHE CR fluxes today, with an anisotropy pattern that reflects the dark 
matter distribution in the halo of our Galaxy. 

It is interesting to note that one of the prime motivations of the inflationary paradigm 
was to dilute excessive production of "dangerous relics" such as topological defects and su- 
perheavy stable particles. However, such objects can be produced right after inflation during 
reheating in cosmologically interesting abundances, and with a mass scale roughly given by 
the inflationary scale which in turn is fixed by the CMB anisotropics to ~ 10 13 GeV [153 



The reader will realize that this mass scale is somewhat above the highest energies observed 
in CRs, which implies that the decay products of these primordial relics could well have 
something to do with EHECRs which in turn can probe such scenarios! 

For dimensional reasons the spatially averaged X particle injection rate can only depend 
on the mass scale mx and on cosmic time t in the combination 

h x (t) =Km p x t- i+p , (13) 

where k and p are dimensionless constants whose value depend on the specific top-down 



scenario 



73], For example, the case p = 1 is representative of scenarios involving release of 



X particles from topological defects, such as ordinary cosmic strings |154 |, necklaces ||155| 



and magnetic monopoles ||156||. This can be easily seen as follows: The energy density p, 



■V 



in a network of defects has to scale roughly as the critical density, p s oc p cr i t oc t 2 , where 



16 



t is cosmic time, otherwise the defects would either start to overdose the Universe, or end 
up having a negligible contribution to the total energy density. In order to maintain this 
scaling, the defect network has to release energy with a rate given by p s = —ap s /t oc t~ 3 , 
where a = 1 in the radiation dominated aera, and a = 2/3 during matter domination. If 
most of this energy goes into emission of X particles, then typically k ~ 0(1). In the flux 
calculations for TD models presented in this paper, it was assumed that the X particles are 
nonrelativistic at decay. 

The X particles could be gauge bosons, Higgs bosons, superheavy fermions, etc. depending 
on the specific GUT. They would have a mass mx comparable to the symmetry breaking 
scale and would decay into leptons and/or quarks of roughly comparable energy. The quarks 
interact strongly and hadronize into nucleons (Ns) and pions, the latter decaying in turn 
into 7-rays, electrons, and neutrinos. Given the X particle production rate, dnx/dt, the 
effective injection spectrum of particle species a (a = , y,N,e ± ,i>) via the hadronic channel 
can be written as (dnx / dt)(2/mx)(dN a / dx) , where x = 2E/mx, and dN a /dx is the relevant 
fragmentation function (FF). 

We adopt the Local Parton Hadron Duality (LPHD) approximation [157] according to 
which the total hadronic FF, dN^/dx, is taken to be proportional to the spectrum of the 
partons (quarks/gluons) in the parton cascade (which is initiated by the quark through 
perturbative QCD processes) after evolving the parton cascade to a stage where the typical 
transverse momentum transfer in the QCD cascading processes has come down to ~ -R -1 ~ 
few hundred MeV, where R is a typical hadron size. The parton spectrum is obtained 
from solutions of the standard QCD evolution equations in modified leading logarithmic 



approximation (MLLA) which provides good fits to accelerator data at LEP energies ||157 



We will specifically use a recently suggested generalization of the MLLA spectrum that 
includes the effects of supersymmetry [ |158|| . Within the LPHD hypothesis, the pions and 
nucleons after hadronization have essentially the same spectrum. The LPHD does not, 
however, fix the relative abundance of pions and nucleons after hadronization. Motivated by 
accelerator data, we assume the nucleon content /jv of the hadrons to be in the range 3 to 
10%, independent of energy, and the rest pions distributed equally among the three charge 
states. For x < 0.1 this ratio as well as the spectra are roughly consistent with recent more 



detailed Monte Carlo simulations ||159| , |160|| . This is also the range relevant for comparison 
with observations as long as mx 3> 10 12 GeV, which covers the range considered in the 
present paper. Resulting normalization ambiguities ||159|| are of factors of a few and thus 



subdominant in light of many other uncertainties at the current stage of these scenarios. 
The standard pion decay spectra then give the injection spectra of 7-rays, electrons, and 
neutrinos. For more details concerning uncertainties in the X particle decay spectra see 



Ref. [82 



4.2 Numerical Simulations 

The 7-rays and electrons produced by X particle decay initiate electromagnetic (EM) cas- 
cades on low energy radiation fields such as the CMB. The high energy photons undergo 
electron-positron pair production (PP; 77^ — > e~e + ), and at energies below ~ 10 14 eV they 
interact mainly with the universal infrared and optical (IR/O) backgrounds, while above 
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~ 100 EeV they interact mainly with the universal radio background (URB). In the Klein- 
Nishina regime, where the CM energy is large compared to the electron mass, one of the out- 
going particles usually carries most of the initial energy. This "leading" electron (positron) 
in turn can transfer almost all of its energy to a background photon via inverse Compton 
scattering (ICS; cjb —>■ e'7). EM cascades are driven by this cycle of PP and ICS. The 
energy degradation of the "leading" particle in this cycle is slow, whereas the total number 
of particles grows exponentially with time. This makes a standard Monte Carlo treatment 
difficult. Implicit numerical schemes have therefore been used to solve the relevant kinetic 
equations. A detailed account of the transport equation approach used in the calculations 



whose results are presented in this contribution can be found in Ref. ||161|| . All EM interac- 
tions that influence the 7-ray spectrum in the energy range 10 s eV < E < 10 25 eV, namely 
PP, ICS, triplet pair production (TPP; t^ h — > ee~e + ), and double pair production (DPP, 
776 - * e~e + e~e + ), as well as synchrotron losses of electrons in the large scale extragalactic 
magnetic field (EGMF), are included. 

Similarly to photons, UHE neutrinos give rise to neutrino cascades in the primordial 
neutrino background via exchange of W and Z bosons ||162| , |7Tf . Besides the secondary neu- 
trinos which drive the neutrino cascade, the W and Z decay products include charged leptons 
and quarks which in turn feed into the EM and hadronic channels. Neutrino interactions 
become especially significant if the relic neutrinos have masses m v in the eV range and thus 
constitute hot dark matter, because the Z boson resonance then occurs at an UHE neutrino 
energy E res = 4 x 10 21 (eV/m„) eV. In fact, this has been proposed as a significant source of 
EHECRs [78, 80]. Motivated by recent experimental evidence for neutrino mass we assume 



the neutrino masses to be m Ve = 0.1 eV, = m Vr = 1 eV and implemented the rele- 
vant W boson interactions in the t-channel and the Z boson exchange via t- and s-channel. 
Hot dark matter is also expected to cluster, potentially increasing secondary 7-ray and nu- 
cleon production |78, [80]. This influences mostly scenarios where X decays into neutrinos 



only. We parametrize massive neutrino clustering by a length scale l v and an overdensity 
/„ over the average density n u . The Fermi distribution with a velocity dispersion v yields 
fu ^ v 3 ml/(2ir) 3 / 2 /n u ~ 330 (v/500 kmsec" 1 ) 3 {m v /eVf Therefore, values of l v ~ 

few Mpc and f v ~ 20 are conceivable on the local Supercluster scale ^fj . 

The relevant nucleon interactions implemented are pair production by protons (p'-fb - > 
pe~e + ), photoproduction of single or multiple pions (N^b — > N rnr, n > 1), and neutron 
decay. In TD scenarios, the particle injection spectrum is generally dominated by the "pri- 
mary" 7-rays and neutrinos over nucleons. These primary 7-rays and neutrinos are produced 
by the decay of the primary pions resulting from the hadronization of quarks that come from 
the decay of the X particles. The contribution of secondary 7-rays, electrons, and neutrinos 
from decaying pions that are subsequently produced by the interactions of nucleons with 
the CMB, is in general negligible compared to that of the primary particles; we nevertheless 
include the contribution of the secondary particles in our code. 

In principle, new interactions such as the ones involving a TeV quantum gravity scale 
can not only modify the interactions of primary particles in the detector, as discussed in 
Sect. 3.3, but also their propagation. However, 7— rays and nuclei interact mostly with 
the CMB and IR for which the CM energy is at most ~ 30(£/10 14 GeV) 1/2 GeV. At such 
energies the new interactions are much weaker than the dominating electromagnetic and 
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strong interactions. It has been suggested recently ||164|| that new interactions may notably 



influence UHE neutrino propagation. However, for neutrinos of mass m v , the CM energy is 
~ 100(£'/10 14 GeV) 1//2 (m l/ /0.1 eV) 1//2 GeV, and therefore, UHE neutrino propagation would 
only be significantly modified for neutrino masses significantly larger than 0.1 eV. We will 
therefore ignore this possibility here. 

We assume a flat Universe with no cosmological constant, and a Hubble constant of 
h = 0.65 in units of 100 km sec _1 Mpc _1 throughout. The numerical calculations follow 
all produced particles in the EM, hadronic, and neutrino channel, whereas the often-used 
continuous energy loss (CEL) approximation (e.g., ||165|| ) follows only the leading cascade 



particles. The CEL approximation can significantly underestimate the cascade flux at lower 
energies. 

The two major uncertainties in the particle transport are the intensity and spectrum of 
the URB for which there exists only an estimate above a few MHz frequency ||166|| , and the 



average value of the EGMF. To bracket these uncertainties, simulations have been performed 
for the observational URB estimate from Ref. | |166| | that has a low-frequency cutoff at 2 



MHz ("minimal"), and the medium and maximal theoretical estimates from Ref. |83| , as 
well as for EGMFs between zero and 10 -9 G, the latter motivated by limits from Faraday 



rotation measurements ||167|| . A strong URB tends to suppress the UHE 7-ray flux by 
direct absorption whereas a strong EGMF blocks EM cascading (which otherwise develops 
efficiently especially in a low URB) by synchrotron cooling of the electrons. For the IR/O 
background we used the most recent data [ |168|| . 

4.3 Results: 7— ray and Nucleon Fluxes 

Fig. |H shows results from Ref. |82[ for the time averaged 7— ray and nucleon fluxes in a 
typical TD scenario, assuming no EGMF, along with current observational constraints on 
the 7— ray flux. The spectrum was optimally normalized to allow for an explanation of the 
observed EHECR events, assuming their consistency with a nucleon or 7— ray primary. The 
flux below < 2 x 10 19 eV is presumably due to conventional acceleration in astrophysical 
sources and was not fit. Similar spectral shapes have been obtained in Ref. [ 169j| , where 



the normalization was chosen to match the observed differential flux at 3 x 10 20 eV. This 
normalization, however, leads to an overproduction of the integral flux at higher energies, 
whereas above 10 20 eV, the fits shown in Figs. ^| and [7| have likelihood significances above 
50% (see Ref. ||170|| for details) and are consistent with the integral flux above 3 x 10 20 eV 



estimated in Refs. [T2|, [13]. The PP process on the CMB depletes the photon flux above 
100 TeV, and the same process on the IR/O background causes depletion of the photon flux 
in the range 100 GeV-100 TeV, recycling the absorbed energies to energies below 100 GeV 
through EM cascading (see Fig. |). The predicted background is not very sensitive to the 
specific IR/O background model, however ||171| . The scenario in Fig. || obviously obeys all 



current constraints within the normalization ambiguities and is therefore quite viable. Note 
that the diffuse 7— ray background measured by EGRET up to 10 GeV puts a strong 
constraint on these scenarios, especially if there is already a significant contribution to this 



background from conventional sources such as unresolved 7— ray blazars [172]. However, the 



7— ray background constraint can be circumvented by assuming that TDs or the decaying 
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long lived X particles do not have a uniform density throughout the Universe but cluster 
within galaxies ||173|| . As can also be seen, at energies above 100 GeV, TD models are not 
significantly constrained by observed 7— ray fluxes yet (see Ref. |19| for more details on these 
measurements). 

Fig. [7] shows results for the same TD scenario as in Fig. |6], but for a high EGMF ~ 10~ 9 G, 
somewhat below the current upper limit | |167| |. In this case, rapid synchrotron cooling of 
the initial cascade pairs quickly transfers energy out of the UHE range. The UHE 7— ray 
flux then depends mainly on the absorption length due to pair production and is typically 
much lower ||165| , |174|| . (Note, though, that for mx ^ 10 25 eV, the synchrotron radiation 
from these pairs can be above 10 20 eV, and the UHE flux is then not as low as one might 
expect.) We note, however, that the constraints from the EGRET measurements do not 
change significantly with the EGMF strength as long as the nucleon flux is comparable 
to the 7— ray flux at the highest energies, as is the case in Figs. | and |7|. The results of 



Ref. [82 differ from those of Ref. |169| which obtained more stringent constraints on TD 



models because of the use of an older fragmentation function from Ref. ||175|| , and a stronger 



dependence on the EGMF because of the use of a weaker EGMF which lead to a dominance 
of 7— rays above ~ 10 20 eV. 

The energy loss and absorption lengths for UHE nucleons and photons are short (< 100 
Mpc). Thus, their predicted UHE fluxes are independent of cosmological evolution. The 
7— ray flux below ~ 10 11 eV, however, scales as the total X particle energy release integrated 
over all redshifts and increases with decreasing p | 176 |. For mx = 2 x 10 16 GeV, scenarios 
with p < 1 are therefore ruled out (as can be inferred from Figs. |6| and ^), whereas constant 
comoving injection models (p = 2) are well within the limits. 

We now turn to signatures of TD models at UHE. The full cascade calculations predict 
7— ray fluxes below 100 EeV that are a factor ~ 3 and ~ 10 higher than those obtained 
using the CEL or absorption approximation often used in the literature, in the case of strong 
and weak URB, respectively. Again, this shows the importance of non- leading particles in 
the development of unsaturated EM cascades at energies below ~ 10 22 eV. Our numerical 
simulations give a 7/CR flux ratio at 10 19 eV of ~ 0.1. The experimental exposure required 
to detect a 7— ray flux at that level is ~ 4 x 10 19 cm 2 sec sr, about a factor 10 smaller than 
the current total experimental exposure. These exposures are well within reach of the Pierre 



Auger Cosmic Ray Observatories |Bfl , which may be able to detect a neutral CR component 
down to a level of 1% of the total flux. In contrast, if the EGMF exceeds ~ 10~ n G, then 
UHE cascading is inhibited, resulting in a lower UHE 7— ray spectrum. In the 10~ 9 G 
scenario of Fig. [7|, the 7/CR flux ratio at 10 19 eV is 0.02, significantly lower than for no 
EGMF. 

It is clear from the above discussions that the predicted particle fluxes in the TD scenario 
are currently uncertain to a large extent due to particle physics uncertainties (e.g., mass and 
decay modes of the X particles, the quark fragmentation function, the nucleon fraction fx, 
and so on) as well as astrophysical uncertainties (e.g., strengths of the radio and infrared 
backgrounds, extragalactic magnetic fields, etc.). More details on the dependence of the 
predicted UHE particle spectra and composition on these particle physics and astrophysical 
uncertainties are contained in Ref. |32 . A detailed study of the uncertainties involved in the 



propagation of UHE nucleons, 7— rays, and neutrinos is currently underway [ |177| |. 
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We stress here that there are viable TD scenarios which predict nucleoli fluxes that 
are comparable to or even higher than the 7— ray flux at all energies, even though 7— rays 
dominate at production. This occurs, e.g., in the case of high URB and/or for a strong 
EGMF, and a nucleon fragmentation fraction of ~ 10%; see, for example, Fig. [7| Some 
of these TD scenarios would therefore remain viable even if EHECR induced EAS should 
be proven inconsistent with photon primaries (see, e.g., Ref. | |1 73 ] ) . This is in contrast 
to scenarios with decaying massive dark matter in the Galactic halo which, due to the 
lack of absorption, predict compositions directly given by the fragmentation function, i.e. 
domination by 7— rays. 

The normalization procedure to the EHECR flux described above imposes the constraint 
Qeuecr ~ 10~ 22 eV cm -3 sec -1 within a factor of a few ||169| , p2| , |179|| for the total energy 



release rate Qo from TDs at the current epoch. In most TD models, because of the unknown 
values of the parameters involved, it is currently not possible to calculate the exact value 
of Qo from first principles, although it has been shown that the required values of Qo (in 
order to explain the EHECR flux) mentioned above are quite possible for certain kinds of 
TDs. Some cosmic string simulations and the necklace scenario suggest that defects may 
lose most of their energy in the form of X particles and estimates of this rate have been 
given ||180| , |155|| . If that is the case, the constraint on Qehecr translates via Eq. (|T3| ) into 



a limit on the symmetry breaking scale r\ and hence on the mass mx of the X particle: 
7] ~ mx ^ 10 13 GeV ||181|| . Independently of whether or not this scenario explains EHECR, 
the EGRET measurement of the diffuse GeV 7— ray background leads to a similar bound, 
Qem ~ 2.2 x 1CT 23 h(3p — 1) eVcm -3 sec -1 , which leaves the bound on rj and mx practically 
unchanged. Furthermore, constraints from limits on CMB distortions and light element 
abundances from 4 He-photodisintegration are comparable to the bound from the directly 
observed diffuse GeV 7-rays ||1 76|| . That these crude normalizations lead to values of rj in the 
right range suggests that defect models require less fine tuning than decay rates in scenarios 
of metastable massive dark matter. 



4.4 Results: Neutrino Fluxes 

As discussed in Sect. 4.1, in TD scenarios most of the energy is released in the form of 
EM particles and neutrinos. If the X particles decay into a quark and a lepton, the quark 
hadronizes mostly into pions and the ratio of energy release into the neutrino versus EM 
channel is r ~ 0.3. 

Fig. [5] shows predictions of the total neutrino flux for the same TD model on which Fig. |6] 
is based. In the absence of neutrino oscillations the electron neutrino and anti-neutrino fluxes 
are about a factor of 2 smaller than the muon neutrino and anti-neutrino fluxes, whereas the 
t— neutrino flux is in general negligible. In contrast, if the interpretation of the atmospheric 
neutrino deficit in terms of nearly maximal mixing of muon and r— neutrinos proves correct, 
the muon neutrino fluxes shown in Fig. [| would be maximally mixed with the r— neutrino 
fluxes. The TD flux component clearly dominates above ~ 10 19 eV. 

In order to translate neutrino fluxes into event rates, one has to fold in the interaction 
cross sections with matter. At UHEs these cross sections are not directly accessible to 
laboratory measurements. Resulting uncertainties therefore translate directly to bounds on 
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neutrino fluxes derived from, for example, the non-detection of UHE muons produced in 
charged- current interactions. In the following, we will assume the estimate Eq. (0) based on 
the Standard Model for the charged- current muon-neutrino-nucleon cross section <j u n if not 
indicated otherwise. 

For an (energy dependent) ice or water equivalent acceptance A(E) (in units of volume 
times solid angle), one can obtain an approximate expected rate of UHE muons produced 
by neutrinos with energy > E, R(E), by multiplying A(E)cr uN (E)nii 2 o (where n H2 o is the 
nucleon density in water) with the integral muon neutrino flux ~ Ej v . This can be used to 
derive upper limits on diffuse neutrino fluxes from a non-detection of muon induced events. 
Fig. shows bounds obtained from several experiments: The Frejus experiment derived 
upper bounds for E > 10 12 eV from their non-detection of almost horizontal muons with 
an energy loss inside the detector of more than 140 MeV per radiation length [|142|| . The 
AMANDA neutrino telescope has established an upper limit in the TeV-PeV range J52" 



The Fly's Eye experiment derived upper bounds for the energy range between ~ 10 eV 
and ~ 10 20 eV |3B[ from the non-observation of deeply penetrating particles. The NASA 



Goldstone radio telescope has put an upper limit from the non-observation of pulsed radio 
emission from cascades induced by neutrinos above ~ 10 20 eV in the lunar regolith. The 
AKENO group has published an upper bound on the rate of near-horizontal, muon-poor 
air showers ||182|| (not shown in Fig. |5]). Horizontal air showers created by electrons, muons 
or tau leptons that are in turn produced by charged-current reactions of electron, muon or 
tau neutrinos within the atmosphere have recently also been pointed out as an important 



method to constrain or measure UHE neutrino fluxes [53| with next generation detectors. 

Clearly, the SLBY model is not only consistent with the constraints discussed in Sect. 4.3, 
but also with all existing neutrino flux limits within 2-3 orders of magnitude. What, then, 
are the prospects of detecting UHE neutrino fluxes predicted by TD models? In a 1 km 3 2n sr 
size detector, the SLBY scenario from Fig. || for example, predicts a muon-neutrino event 
rate of ~ 0.08 yr _1 , and an electron neutrino event rate of ~ 0.05 yr -1 above 10 19 eV, where 
"backgrounds" from conventional sources should be negligible. Further, the muon-neutrino 
event rate above 1 PeV should be ~ 0.6 yr -1 , which could be interesting if conventional 
sources produce neutrinos at a much smaller flux level. Moreover, the neutrino flux around 
10 17 eV could have a slight enhancement due to neutrinos from muons produced by interac- 
tions of UHE photons and electrons with the CMB at high redshift ||183| , an effect that has 
not been taken into account in the figures shown here. Of course, above ~ 100 TeV, instru- 
ments using ice or water as detector medium, have to look at downward going muon and 
electron events due to neutrino absorption in the Earth. However, r— neutrinos obliterate 
this Earth shadowing effect due to their regeneration from r decays ||106|| . The presence of 
t— neutrinos, for example, due to mixing with muon neutrinos, as suggested by recent ex- 
perimental results from Super-Kamiokande, can therefore lead to an increased upward going 



event rate ||104 |. As mentioned in Sect. 2, r neutrinos skimming the Earth at small angles 



below the horizon can also lead to an increase of sensitivity of fluorescence and ground array 



detectors 64, 65, 66 



For detectors based on the fluorescence technique such as the HiRes j|2J and the Telescope 
Array |43| (see Sect. 2), the sensitivity to UHE neutrinos is often expressed in terms of 
an effective aperture a(E) which is related to A(E) by a(E) = A(E)<j uN (E)ri}i 2 o. F° r 
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the cross section of Eq. (|7|), the apertures given in Ref. |42| for the HiRes correspond to 
A(E) ~ 3 km 3 x 2vrsr for E > 10 19 eV for muon neutrinos. The expected acceptance of the 
ground array component of the Pierre Auger project for horizontal UHE neutrino induced 
events is A(10 19 eV) ~ 20km 3 sr and A(10 23 eV) ~ 200km 3 sr |k§, with a duty cycle close 
to 100%. We conclude that detection of neutrino fluxes predicted by scenarios such as the 
SLBY scenario shown in Fig. |5| requires running a detector of acceptance > 10 km 3 x 27rsr 
over a period of a few years. Apart from optical detection in air, water, or ice, other methods 
such as acoustical and radio detection |32] (see, e.g., the RICE project |58| for the latter) or 



even detection from space |46|, [47|, [48], |4£| appear to be interesting possibilities for detection 
concepts operating at such scales (see Sect. 2). For example, the space based OWL/ Air Watch 
satellite concept would have an aperture of ~ 3 x 10 6 km 2 sr in the atmosphere, corresponding 
to A(E) ~ 6 x 10 4 km 3 sr for E > 10 20 eV, with a duty cycle of ~ 0.08 |§, g|. The 
backgrounds seem to be in general negligible |7l], 184|| . As indicated by the numbers above 



and by the projected sensitivities shown in Fig. [5], the Pierre Auger Project and especially 
the space based AirWatch type projects should be capable of detecting typical TD neutrino 
fluxes. This applies to any detector of acceptance > 100 km 3 sr. Furthermore, a 100 day 
search with a radio telescope of the NASA Goldstone type for pulsed radio emission from 
cascades induced by neutrinos or cosmic rays in the lunar regolith could reach a sensitivity 
comparable or better to the Pierre Auger sensitivity above ~ 10 19 eV [|60|1 . 

A more model independent estimate ||179|| for the average event rate R{E) can be made 



if the underlying scenario is consistent with observational nucleon and 7— ray fluxes and the 
bulk of the energy is released above the PP threshold on the CMB. Let us assume that the 
ratio of energy injected into the neutrino versus EM channel is a constant r. As discussed in 
Sect. 4.3, cascading effectively reprocesses most of the injected EM energy into low energy 
photons whose spectrum peaks at ~ lOGeV ||171|| . Since the ratio r remains roughly un- 
changed during propagation, the height of the corresponding peak in the neutrino spectrum 
should roughly be r times the height of the low-energy 7— ray peak, i.e., we have the condition 



max E 



E ■ j v (E) ~ r maxg [E j 7 (E)] . Imposing the observational upper limit on the diffuse 

7— ray flux around 10 GeV shown in Fig. [|, maxg E 2 j Uf ^(E) < 2 x 10 3 reVcm~ 2 ; 
then bounds the average diffuse neutrino rate above PP threshold on the CMB, giving 



_1 sr 



R(E) < 0.34 r 



A(E) 



1 km x 2n sr 



E 



10 19 eV 



-0.6 



yr- 1 (E> 10 i5 eV), (14) 



assuming the Standard Model cross section Eq. (|7|). Comparing this with the flux bounds 
shown in Fig. |5| results in an upper bound on r. For example, the Fly's Eye bound translates 
into r < 20(.E/10 19 eV) ' 1 . We stress again that TD models are not subject to the Waxman 
Bahcall bound because the nucleons produced are considerably less abundant than and are 
not the primaries of produced 7— rays and neutrinos. 

In typical TD models such as the one discussed above where primary neutrinos are pro- 
duced by pion decay, r ~ 0.3. However, in TD scenarios with r ^> 1 neutrino fluxes are 
only limited by the condition that the secondary 7— ray flux produced by neutrino interac- 
tions with the relic neutrino background be below the experimental limits. In this case the 
observed EHECR flux would be produced by the Z-burst mechanism discussed in Sect. 3.1. 
An example for such a scenario is given by X particles exclusively decaying into neutrinos 
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(although this is not very likely in most particle physics models, but see Ref. |82[] and Fig. || 
for a scenario involving topological defects and Ref. |81] for a scenario involving decaying 
superheavy relic particles, both of which explain the observed EHECR events as secondaries 
of neutrinos interacting with the primordial neutrino background). Such scenarios predict 
appreciable event rates above ~ 10 19 eV in a km 3 scale detector, but require unrealistically 
strong clustering of relic neutrinos (a homogeneous relic neutrino overdensity would make 
the EGRET constraint only more severe because neutrino interactions beyond ~ 50 Mpc 
contribute to the GeV 7— ray background but not to the UHECR flux). A detection would 
thus open the exciting possibility to establish an experimental lower limit on r. Being based 
solely on energy conservation, Eq. (14) holds regardless of whether or not the underlying TD 
mechanism explains the observed EHECR events. 

The transient neutrino event rate could be much higher than Eq. ([14]) in the direction 
to discrete sources which emit particles in bursts. Corresponding pulses in the EHE nucleon 
and 7— ray fluxes would only occur for sources nearer than ~ 100 Mpc and, in case of 
protons, would be delayed and dispersed by deflection in Galactic and extragalactic magnetic 
fields [ 185 , |186|| . The recent observation of a possible clustering of CRs above ~ 4 x 10 19 eV by 
the AG AS A experiment [ |187| ] might suggest sources which burst on a time scale t& <C 1 yr. A 
burst fluence of ~ r A(E) /lkm 3 x 27rsr (E/10 19 eV) -0,6 neutrino induced events within a 
time if, could then be expected. Associated pulses could also be observable in the GeV — TeV 
7— ray flux if the EGMF is smaller than ~ 10 -15 G in a significant fraction of extragalactic 
space 



In contrast to roughly homogeneous sources and/or mechanisms with branching ratios 
r 3> 1, in scenarios involving clustered sources such as metastable superheavy relic particles 
decaying with r ~ 1, the neutrino flux is comparable to (not significantly larger than) 
the UHE photon plus nucleon fluxes and thus comparable to the universal cosmogenic flux 
marked u Nj" in Fig. [^. This can be understood because the neutrino flux is dominated 
by the extragalactic contribution which scales with the extragalactic nucleon and 7— ray 
contribution in exactly the same way as in the unclustered case, whereas the extragalactic 
contribution to the "visible" flux to be normalized to the UHECR data is much smaller 
in the clustered case. The resulting neutrino fluxes in these scenarios would thus be much 
harder to detect even with next generation experiments. 

For recent compilations of UHE neutrino flux predictions from astrophysical and TD 
sources see Refs. [[47], |189|| and references therein. 



5 Conclusions 

Ultra-high energy cosmic rays have the potential to open a window to and act as probes of 
new particle physics beyond the Standard Model as well as processes occurring in the early 
Universe at energies close to the Grand Unification scale. Even if their origin will turn out to 
be attributable to astrophysical shock acceleration with no new physics involved, they will 
still be witnesses of one of the most energetic processes in the Universe. The future appears 
promising and exciting due to the anticipated arrival of several large scale experiments. 
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Figure 4: Fluxes of neutrinos (dashed and dashed-dotted, as indicated), 7— rays (solid), and 
nucleons (dotted) predicted by the Z-burst mechanism for m Ue = 0.1 eV, m v = m VT = leV, 
for homogeneously distributed sources emitting neutrinos with an E~ x spectrum (equal for 
all flavors) up to 10 22 eV and an E~ 2 7— ray spectrum of equal power up to 100 TeV. Injection 
rates were assumed co movingly constant up to z = 2. The relic neutrino overdensity was 
assumed to be 200 over 5 Mpc. The calculation used the code described in Ref. and 



assumed the lower limit of the universal radio background |83| and a vanishing extragalactic 
magnetic field. 1 sigma error bars are the combined data from the Haverah Park ||, the Fly's 
Eye [0 



and the AGASA experiments above 10 19 eV. Also shown are piecewise power 
law fits to the observed charged CR flux (thick solid line) and the EGRET measurement of 
the diffuse 7— ray flux between 30 MeV and 100 GeV |84j] (solid line on left margin). Points 
with arrows represent upper limits on the 7— ray flux from the HEGRA J55| , the Utah- 
Michigan [|S6[] , and the CASA-MIA [^] experiments, and on the neutrino flux from the 
Fly's Eye experiment p8|, the Goldstone radio telescope pOl, and the AMANDA neutrino 



telescope p2"| , as indicated (see text and Ref. [0 for more details). 
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Figure 5: Various neutrino flux predictions and experimental upper limits or projected sen- 
sitivities. Shown are upper limits from the Frejus underground detector ||142|| , the Fly's 
Eye experiment |88[| , the Goldstone radio telescope fl60"f , and the Antarctic Muon and Neu- 
trino Detector Array (AMANDA) neutrino telescope as well as projected neutrino 
flux sensitivities of ICECUBE, the planned kilometer scale extension of AMANDA |[56| , 
the Pierre Auger Project (for electron and tau neutrinos separately) and the proposed 
space based OWL |46| concept. Neutrino fluxes are shown for the atmospheric neutrino 
background [|143|1 (hatched region marked "atmospheric"), for EHECR interactions with the 
CMB ||144|| ( "JV7" , dashed range indicating typical uncertainties for moderate source evolu- 
tion), and for the "top-down" model (marked "SLBY") corresponding to Fig. |j] of Sect. 4 
below, where EHECR and neutrinos are produced by decay of superheavy relics. The top- 
down fluxes are shown for electron-, muon, and tau-neutrinos separately, assuming no (lower 
z^-curve) and maximal — v T mixing (upper ^ T -curve, which would then equal the in- 
flux), respectively. The Waxman-Bahcall bound in the version of Mannheim, Protheroe, 
and Rachen ( "WB/MPR-bound" ) for sources optically thin for the proton primaries, 
and the 7— ray bound ("7— bound") are also shown. 
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Figure 6: Predictions for the differential fluxes of 7— rays (solid line) and protons and neu- 
trons (dotted line) in a TD model characterized by p = 1, m x = 10 16 GeV, and the decay 
mode X — ■> q + q, assuming the supersymmetric modification of the fragmentation func- 
tion ||158|| , with a fraction of about 10% nucleons. The defects have been assumed to be 
homogeneously distributed. The calculation used the code described in Ref. [^] and as- 
sumed an intermediate URB estimate from Ref. |83| and an EGMF 10~ n G. Remaining 
line key as in Fig. f|. 
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Figure 7: Same as Fig. |, but for an EGMF of 1(T 9 G. 
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Figure 8: Flux predictions for a TD model characterized by p = 1, mx = 10 14 GeV, with 
X particles exclusively decaying into neutrino-antineutrino pairs of all flavors (with equal 
branching ratio), assuming neutrino masses m Ve = 0.1 eV, m v = m UT = leV. For neutrino 
clustering, an overdensity of ~ 50 over a scale of l v ~ 5 Mpc was assumed. The calculation 
assumed an intermediate URB estimate from Ref. [83] and an EGMF 10~ n G. The line 
key is as in Fig. |]. 
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